Scaling analysis of the magnetoresistance in Gai-^Mn^As : Evidence for strong 

fluctuation and interaction effects 
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We compare experimental resistivity data on Gai-^Mn^As films with theoretical calculations 
using a scaling theory for strongly disordered ferromagnets. All characteristic features of the tem- 
perature dependence of the resistivity can be quantitatively understood through this approach as 
originating from the close vicinity of the metal-insulator transition. In particular, we find that the 
magnetic field induced changes in resistance cannot be explained within a mean-field treatment of 
the magnetic state, and that accounting for thermal fluctuations is crucial for a quantitative analy- 
sis. Similarly, while the non-interacting scaling theory is in reasonable agreement with the data, we 
find clear evidence in favor of interaction effects at low temperatures. 

PACS numbers: 72.20.i, 72.80.r, 75.50.Pp, 72.20.My 



Over the past decade, Gai_ x Mn x As has been the most 
extensively studied ferromagnetic semiconductor because 
it can be used in proof-of-concept spintronic devices 
In particular, the carrier-mediated ferromagnetism in 
this material makes it attractive for device integration 
with the technologically mature III-V semiconductors 
such as GaAs, where bandgap engineering allows sys- 
tematic modulation of the carrier density in heterostruc- 
ture devices. Furthermore, we note that an all-electrical 
III-V semiconductor spintronic device has been recently 
demonstrated [2]. It is important within this general 
context to develop a fundamental understanding of the 
interplay between carrier transport and magnetism in 
Gai_2;Mn x As . One of the basic but least understood 
properties of Gai_ x Mn x As is the temperature and mag- 
netic field dependence of its resistivity Q • 

For typical Gai_ x Mn x As samples the resistivity in- 
creases with decreasing temperature above the Curie 
temperature, Tc, but then it suddenly drops below Tc 
thereby resulting in a resistivity peak at Tc, which then 
gradually broadens and gets shifted to higher tempera- 
tures upon application of a magnetic field (see Fig. []}. 
In "high quality samples", the aforementioned increase 
is less pronounced and one observes a broad shoulder 
rather then a resistivity peak. The resistivity also shows 
another upturn at much lower temperatures. There have 
been a number of attempts to explain the resistivity peak 
so far 3, which invoked, e.g., scattering by critical fluc- 
tuations 0], the formation of magnetic polarons [!, @], 
'dynamical' mean field calculations [2j, or the interplay 
with universal conductance fluctuations [§]. Although 
these theoretical approaches have been quite successful 
in addressing a particular range or qualitative aspect of 
the data, however, a theoretical framework that could 
quantitatively explain all characteristic features observed 



in Gai_ x Mn x As has not been available so far. Besides 
the simplicity of our model and the clear set of reason- 
able parameter choices we make, the attractive feature of 
the approach presented in this paper stems from the fact 
that it is capable of quantitatively account for all distinct 
features of the experimental data: (i) the gradually in- 
creasing resistance as the temperature is lowered towards 
Tc', (ii) the pronounced peak precisely at Tc; (Hi) the 
upturn in resistance at low temperatures, together with 
the finite resistance intercept for metallic samples; (iv) 
the precise amount with which an external magnetic field 
depletes the resistance peak at Tc and shifts it towards 
higher temperatures; (v) the remarkable "non-crossing" 
constraint of the experimental data, a constraint that 
some other theories fail: the resistance curves that be- 
long to different values of the external field do not cross. 

In this Letter, we present a theoretical analysis of the 
temperature dependent resistivity of Gai_ x Mn x As and 
show good agreement with a range of measured samples 
despite of the simple assumptions we make. First of all, 
we remark that the features mentioned above are reminis- 
cent of localization effects interplaying with magnetism. 
In fact, most /// — V ferromagnetic semiconductors are 
very bad conductors because the charged dopants (Mn for 
Gai-^Mn^As ) introduce very large disorder Q and can 
lead even to the formation of an impurity band 0. As 
we show later, even for annealed and relatively highly- 
doped samples, a naive estimate gives fcpZ ~ 0.3, with 
fcp the Fermi momentum and / the mean free path. This 
value clearly hints towards the importance of disorder, 
and suggests that it is necessary to beyond the pop- 
ular free electron picture frequently used in the litera- 
[T3 |. Strikingly similar resistivity anomalies 
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ture 

have been also observed in other types of magnetic semi- 
conductors [IH , as well as some manganites 14 , 1|| , with 
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various semi-phenomenological frameworks available for 
explaining these phenomena in terms of localization the- 
ory HI, la, 12 1- However, these approaches focus explic- 
itly on the localized phase. Such a starting point is not 
suitable for a fairly large class of experimentally mea- 
sured Gai-^Mn^As samples that are not insulators, but 
very poor metals, close to the localization transition. 

In Ref. [l8| , we developed a scaling theory of magneto- 
resistance, wherein all the characteristic features of the 
resistivity anomalies appear naturally. This previous 
work primarily focused on the localized phase, where 
Mott's variable range hopping formula can be used to 
compute the resistivity. The purpose of the present pa- 
per is to extend this theory to metallic samples, and 
show that it qualitatively and quantitatively explains de- 
tailed aspects of the magneto-transport properties. To 
understand the properties of Gai_ x Mn x As we must take 
into account electron-electron interaction too. Never- 
theless, as explained later, the concept of universal one- 
parameter scaling carries over to this material due to the 
large intrinsic spin-orbit coupling of Gai-^Mn^As and 
the presence of magnetic moments. 

To test the theoretical approach, we have measured a 
series of annealed and unannealed Gai-^Mn^As samples 
grown by molecular beam epitaxy in which the Mn con- 
centration (x) was systematically varied between 0.0135 
and 0.067. The temperature dependent resistivity data 
shown in Fig. [T] were measured in different magnetic fields 
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FIG. 1: Comparison between the experimental data and the 
theoretical results at magnetic fields H =0, 3, 6 and 9 T. 
Dots represent experimental data, solid lines are theoretical 
fits. Each figure indicates the value of q' that is used for 
theoretical fitting. 



inside a commercial cryostat (Quantum Design PPMS), 
with the magnetic field normal to the sample plane. De- 
tails of the sample growth, materials characterization and 
annealing protocol were reported elsewhere [ll| • 

In the present paper we shall use a slightly modified 
version of the scaling approach applied for disordered 



conductors. [2fJ, |21|, |22| In Ref. [18| we considered non- 
interacting electrons, where a scaling theory can be con- 
structed in terms of the dimensionless conductance, g 
and a lengthscale at which electrons loose their coher- 
ence. In the presence of interactions, one also needs to 
introduce the dimensionless interaction parameters in the 
triplet and singlet channels, jt and j s , respectively. How- 
ever, for Gai-^Mn^As an important simplification oc- 
curs: Gai-zMn^As has a very large intrinsic spin-orbit 
gap, A so ~ 4000K. Furthermore, in the ferromagnetic 
phase and in the vicinity of the Curie temperature, Tq, 
the almost classical 5 = 5/2 spins of the M n ions fluctu- 
ate slowly in time, and therefore, at the time scales and 
temperatures of interest, time reversal is also locally bro- 
ken, even in the paramagnetic phase. As a result, as we 
showed already in Ref. |18|, Gai-^Mn^As belongs to the 
unitary class .[2 3] Then the 7t plays no role, and 7 S can 

As a result, the scaling of 



also be set to 7 S = l.[21 



the dimensionless conductance is described by a scaling 
equation, 



d In g 

dx 



/?(<?), 



(1) 



where x = m(£) is a scaling variable with £ = £,(T) a 
lengthscale, and due to the simplifications above, the j3- 
function depends only on g itself [2l|, [24| . In three dimen- 
sions, there is a metal-insulator transition characterized 
by P (9c) = 0, with gc the critical conductance. For 
g > gc and (3 (gc) > 0, the conductor is metallic, and 
the dimensionless conductance increases with increasing 
system size, while (3 (g < gc) < 0, and one recovers an 
insulating state. 

Would we know the /^-function, we could compute the 
resistivity as follows: Suppose we know the typical di- 
mensionless conductance go at an energy scale Tq and at 
the corresponding microscopic length scale, £o = £(7o). 
Then the resistivity of a large three-dimensional conduc- 
tor can be computed by integrating the scaling equation 
up to a length scale x — ln(^(T)/^ ) and cutting the 
system to small cubes of size £ = £ (T) to give 



Q(T,H) = 



Z(t) gc 



' 9(i(T),g /gc) 



(2) 



Note that here we compute the typical conductance. Cor- 
rections due to universal conductance fluctuations can 
give a singular contribution, Q. However, these correc- 
tions are very small for magnets with a short mean free 
path such as Gai-^Mn^As , and therefore they cannot 
explain the resistivity anomaly in Gai_a;Mn x As Q. 
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To compute we need to know the beta-function. 
Unfortunately, while this can and has been determined 
numerically for a non-interacting unitary system [Lsj , it is 
not known for interacting electrons. Here we shall use its 
asymptotic form on the metallic side, /3(g) = 1 — 9q /g, 
with the critical value of the conductance to lowest 
order in the epsilon expansion. 

Furthermore, we need to know the connection between 
the scale £(T) and the temperature T. This connection 
can be established by looking at the pole structure of the 
diffusion propagator, £ 2 (T) - D(T)/(T z(T)), and using 
the Einstein relation, a(T) ~ {e 2 /h) N(0)D(T) together 



with Eq. ([2]) 2l|, |22|. However, since we are interested in 
the metallic regime, with a good approximation, we can 
neglect the energy-dependence of the factor z(T) and set 
it to z = 1. As a result, we have AT(0) £ 3 ~ g(£)/T [Hj], 
which can equally be written as 

5(£/£oj So) T 



fjo 



T 



(3) 



with To the energy scale corresponding to the length scale 

With the (3(g) and £(T) at hand, we then only need one 
more ingredient, the microscopic resistivity, go = go{h, t), 
with t = T/Tq and h = g H S/Tq denoting the di- 
mensionless temperature and magnetic field, respectively. 
In the following, we shall assume that go depends exclu- 
sively on the magnetization m of the sample, 

go (t, h) = g {m(t, h)) « g (l + qm 2 (t, h)) (4) 

where go is the conductance of the unpolarized system. 
This approximation is well justified within a mean-field 
description of the scattering on spin disorder [18j, but it 
also emerges quite naturally for mechanisms where the 
primary role of the local magnetic moments is just to 
polarize the charge carriers and indirectly influence their 
scattering rate |25j . In the present formalism, however, 
the precise microscopic origin of the m-dependence of go 
is of secondary importance. 

The quadratic term in Eq. ([4]) is the leading contribu- 
tion allowed by time reversal symmetry, and provides a 
very good approximation even in the extreme case of the 
infinite coupling disordered Kondo lattice [l8| . From the 
fits we found that the parameter q w 0.5 -j- 0.7 for all 
samples, in rough agreement with the results of (25[. 

In principle, we could now use the experimental m(t, h) 
curves to compute the magnetoresistance of a sam- 
ple. Unfortunately, experimentally, at large magnetic 
fields we could not separate the magnetization of the 
Gai-xMn^As film and that of the paramagnetic sub- 
strate, and only magnetization curves at small magnetic 
fields (H = 50 Oe) were available. Therefore, instead of 
using the experimental data, we determined the magne- 
tization m(t, h) in a finite field by performing simulations 
for a diluted spin system. To simulate Gai_ x Mn-jAs we 
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FIG. 2: Temperature dependence of the normalized mag- 
netization in an external (H=50 Oe) in-plane magnetic field. 
For each sample the we rescaled T by the corresponding Tc- 
The solid line is the result of a Monte Carlo calculation for a 
diluted magnetic semiconductor Gai_ I Mn I As by assuming a 
RKKY interaction between the localized Mn spins. 



placed magnetic ions with a given concentration at ran- 
dom positions of an FCC lattice following the procedure 
of Ref. [26|| . We assumed an RKKY interaction between 
the Mn spins and computed the magnetization curves 
m(t, h) by performing a Monte Carlo simulation. This 
procedure reproduces the m(t, h) curves, which turn out 
to be almost concentration-independent (apart from the 
limit of very small Mn concentrations), and fit very nicely 
the experimentally measured magnetization for h = too 
(see Fig. [5]). Although these magnetization curves look 
rather similar to the ones obtained from a simple mean 
field approximation, Tc is suppressed by a factor ~ 2 
compared to the mean field Curie temperature due to 
thermal fluctuations. These thermal fluctuations being 
very sensitive to the magnetic field, the magnetization 
induced by an external field of ~ 1 T at T = Tc was 
about a factor ~ 1.5 larger than the mean field estimate. 
This increased response was important to obtain the cor- 
rect position of the resistivity maximum in finite mag- 
netic fields. Note, however, that the the curves m(t, h) 
obtained this way have no free fitting parameters. 

We now have all ingredients to compute the mag- 
netoresistance. The temperature and magnetic field- 
dependence of the resistivity then originates from the 
temperature and magnetic field-dependence of the mi- 
croscopic conductance go and that of the scale £: The 
correlation length £ becomes larger as T — + 0, and there- 
fore the resistivity increases. This results ultimately in 
the low-temperature upturn of the resistivity and is also 
responsible for the upturn of the resistivity above Tc ■ At 
very low temperatures this results in a ~ VT-dependencc 
[27j ]. Entering the ferromagnetic phase, or polarizing the 
Mn moments with an external field, on the other hand, 
increases go, and hence decreases the resistivity. It is the 
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TABLE I: Characteristic parameters of the samples analyzed: x is the Mn concentration, g(Tc) is the resistivity at To = Tc, 
and fci? stands for the Fermi momentum obtained by assuming a compensation of 50%. We computed k F l from the Drude 
formula. We also show the Fermi wavelength Xf of the non-magnetic system, the fitted correlation length at To = Tc [£o] 
and at IK [I*(1K)], respectively, the dephasing length £ Drude (Tc) obtained using the Drude estimate, as described in the text. 
Finally we list the fitted values of go/gc- 
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competition of these two effects that yields the resistivity 
anomaly at Tc- 

Eqs. © and (g]) together with Eq. Q provide a full 
theoretical description of the magneto-resistance in terms 
of three parameters for every sample, £o, go/gc, and the 
phenomenological parameter q. Fig. [1] shows the best fits 
obtained in this way for six different samples. For all 
samples we defined £o and go as the scale and dimcn- 
sionless conductance at To = Tc- The parameters of the 
samples are summarized in Table |U Note that the posi- 
tion, the shift and the amplitude of both the resistivity 
maxima at Tc as well as that of the low-temperature 
anomaly work out very nicely for these samples. This 
proves indirectly, that both anomalies are related to the 
vicinity of the metal-insulator transition. 

In Table U we also enumerate the values of UfI we ob- 
tain from the resistivity by assuming a valence hole band 
of effective mass m* = 0.45 m e and a compensation of 
50%. Clearly, the values obtained in this way are incon- 
sistent with a weakly disordered free electron picture, but 
could naturally be explained through the presence of an 
impurity band with an enhanced carrier mass.(l2l|. 

Importantly, the fitted values of £ are smaller than the 
thickness of the films, W w 123 nm even at T = IK, 
and they are in good agreement with the values obtained 
in Refs. [3]. These samples are thus three-dimensional 
from the point of view of conductance properties down 
to these temperatures. Also, £ remains larger than the 
typical Mn-Mn separation ~ lnm over the whole range 
of temperatures, thereby justifying the scaling approach 
used here. In Table Q] we also included the theoretical 
estimate of £ Drude (T) = yj Dn ru de/T where we estimated 
-C>Drude by using the Drude formula and the density of 
states of a parabolic valence band with a renormalizcd 
mass. The values obtained in this way do not depend too 
much on the specific sample, and, apart from an overall 
factor, are in rough agreement with the values extracted 
from the fits. This is somewhat surprising in view of the 
extremely small values of kpl. 

In conclusion, we presented here a systematic study 
of the resistivity of various Gai-^Mn^As samples. We 



argued that even the annealed samples are very close to 
the localization transition and showed that the complete 
magnetic field dependence of the resistivity anomaly at 
Tc as well as the low temperature upturn of the resis- 
tivity can be quantitatively described in terms of a per- 
turbative scaling theory of localization, combined with 
Monte Carlo simulations and a reasonable choice of ma- 
terials parameters. 
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